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We study the long range behavior of a gas whose partition function depends on a parameter 
q and it has been claimed to be a good approximation to the partition function proposed in the 
formulation of nonextensive statistical mechanics. We compare our results, at large temperatures 
' and at the critical point, with the case of Boltzmann-Gibbs thermodynamics for the case of a Bose- 

I Einstein gas. In particular, we find that for all temperatures the long range correlations in a Bose 

. gas decrease when the value of q departs from the standard value q = 1. 
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The formulation of nonextensive statistical mechanics has raised numerous questions regarding its relevance to 
systems with long range interactions, long range microscopic memory, or multifractal properties This formalism 
can also be understood as a generalization of Boltzmann-Gibbs statistics, opening the possibility to get a theoretical 
insight on the thermodynamics of systems whose behavior depart from Boltzmann-Gibbs statistical mechanics. The 
generalized entropy 

S, = Z^\^-l^P'n]' (1) 



is a function of the probability pn for the ensemble to be in the state R and a real parameter q. Equation (ill) becomes 



^ : the Shannon entropy as g 1. In addition, Tsallis' entropy shares all the properties of the Shannon entropy except 
•I— I 

' that of additivity. Considering two independent systems S and E', the entropy Sq satisfies the pseudoadditivity 
property 

qEUE' cS cS' cE oE' 

The probability distribution that results from extremizing the entropy with the constraints 

R 

and 
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is given by the equation 



with the partition function 



[l + Piq-l)iEn~f,N)f^'-''^ 

PR = y > (5) 



Z, = Y,[1 + P{q - l)iEn - ^iV)]i/(^-') , (6) 

R 

and the total energy En — njEj. Due to the mathematical complexity of the partition function in Equation (|^), 
it has been necessary to study the validity of certain approximation procedures. Thus, the study of the consequences 
of nonextensivity has been mainly focused on either assuming a value of q w 1 or by approximating the partition 
function in Eq. by a factorized partition function Q 

^ = nE [1 + /3(9 -m)]'/^'-'^ (7) 

The factorization approach has been shown to be a good approximation to Tsallis partition function outside 
certain narrow temperature interval that shifts to higher values of T when the number of energy levels increases. 
An application to the Ising model and black body radiation can be found in Refs. ||^ and [Q respectively. In this 
approximation the average number of particles with energy e is given by the function 

(n) = 771— ' (8) 

where a = 0,— 1,+1 for Maxwell-Boltzmann, Bose-Einstein and Fermi-Dirac cases respectively. It is important 
to remark that the function (n) has also been obtained by extremizing the entropy related to the generalized 
dimensions of a fractal set. It has been also pointed out that nonextensive thermodynamics could also be 
understood in terms of g-deformations and possibly with the theory of quantum groups. Along this line of work, in 
Ref. we made a study of the basic thermodynamics that result from the particle distribution functions in Eq. (|^) 
for classical and quantum gases. In that article we shown that the high temperature behavior is consistent with the 
thermodynamical limit provided that the internal energy is calculated according to the equation 

(C/) = ^^°°|^Wp)>Vdp. (9) 

By introducing the usual creation and annihilation operator formalism we found that the boson hamiltonian that 
leads to the particle distribution, for a = —1, in Eq.(||) is written as 

i: = ^(£,-M)0,</)„ (10) 

where the operator and its adjoint satisfy a deformed boson algebra. However, a comparison of the heat capacity 
and entropy functions, in Ref. for systems with a particle distribution function as given by Eq. (||) with those for 



Bose and Fermi gases described by quantum group invariant hamiltonians shows that nonextensivity is unrelated 
to quantum group invariance. 

In this paper we calculate the correlation function for boson systems with a particle distribution function as given 
by Eq. (||). Our main motivation in studying this system is twofold. First, it is of theoretical interest to study a 
thermodynamics system which obeys a statistical mechanics that generalizes Boltzmann-Gibbs statistics. Second, 
a calculation of the correlation functions will give us an insight on the long-range behavior dependence of these 
thermodynamical systems on the parameter q. Our calculations will show in a concrete fashion the relation between 
q and long-range behavior, and will indicate whether the thermodynamics resulting from Eq. (^), proposed in Ref. 

has the long range behavior expected to be present in nonextensive thermodynamics. In Sec. |^we calculate the 
correlation function for q =/= I, and discuss the results for some particular values of this parameter. We specialize our 



discussion to the behavior near the critical temperature Tc and at large values of T. In Sec. Ill we summarize our 
results. 

II. CORRELATION FUNCTIONS 

According to our previous work the parameter q cannot be any real number but its values are restricted to 
those such that l/(g— 1) is an integer. In addition, the thermodynamic functions are well defined only in the interval 
1 < g < 1.5. The correlation for a Bose gas with particle distribution according to Eq.(^) is given by 

""^""^^^^[l + Piq-lKe-.r^^'^^-l ^''^ 
As usual, this summation can be approximated by an integration over k, leading after integration over the angles for 
d = 3 to the equation 

where the new variable — (3{q— l)k^h^ /2m and, as usual, a = — and A is the thermal wavelength. The integral 
in Eq. (H) is tabulated |l|| 

where —1 < Re v < i?e(2/i -I- 3/2), c > and & > 0. It is simple to check that all these inequalities are satisfied for 
all the allowed values of q in the interval 1 < 9 < 1.5. Defining a length parameter x < 



^ = ^V i+l('/^i) ' 



and by use of the integral representation 



the correlation function reduces to the expression 



1 



{R/x? 



n/{q-l)-l 



dt. 



(16) 



Xl + a{q^l)){t^ + {R/xY)_ 

Once we perform the summation we get a set of integrals in terms of elementary functions. These integrals are all 
tabulated 0| and the correlation function for some values of q follows 



G(R) = < 



(l/i?A2) [e-«/« - e-V(«/C)'+32-(i?/A)^ _ 2e-^' sin(B')j for q = 5/4 
{l/R\^)[e-^/^ -e-'^cos{B)-^/ie-'^s\n{B)\ for q = 4/3 , 

for q = 3/2 



(l/i?A2) er^/^ - 



(17) 



where ^ = X/2^'Ka is the correlation length and the exponents 



= 127r(i?/A)2 [(l/2)v/(3/2 + a/3)2 + (3/4) + (3/2 + a/3) 
- 127r(i?/A)2 [(1/2)^(3/2 + a/3)2 + (3/4)- (3/2 + a/3) 



= 167r(i?/A)2 (l/2)V(l + a/4)2 + l + (l + a/4) 



B'2 = 167r(i?/A)2 (l/2)v/(l + a/4)2 + l-(l + a/4) 



From these equations it is clear that the correlation functions for q ^ 1 have a leading term 



G(R) « (l/i?A^)e 



minus some smaller additional terms which are absent in the Bose-Einstein case. 



A. Correlation length 



(18) 



In order to compare these correlation functions with the standard correlation function Gs£;(R) wc first look at the 
case of high temperature T >> Tc- Although Eq. (ITsI) has the same functional relation than 



the parameter a is the function 



9-1 



with 



GbeCR) ~ {l/RX')e 



0F(n)(g- 1)3/2A' 



:S2{q) 



2(9-l)/(5-3g)- 



1 °° (_T\m 



1 



3/2 - m - l/{q-l)' 



(19) 



(20) 



(21) 
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In terms of the critical temperature Tc the correlation length ^ for high temperatures is given by the expression 

1/2 



e V 9-1 VG3/2(l,g); \T, 



1/3 / rp \ 1/2 



1 (n) \''WT 

= 2 



J, ^3/2 X 2(9-l)/(5-3g)' 



(22) 



where the function 03/2(2; = 1, (?) was calculated in Ref. In particular 03/2(1, 1) = 2.612 and 63/2(1, q) increases 
with the value of q. A simple inspection shows that the correlation length ^ is smaller than the correlation length for 
the standard Bose-Einstein case 

' 2V^ ( M\ ( _L\ \2M2{T/T!?^f/^] , (23) 



V \^ 2.612 ; \TBE ^ 

where < T^^. Fi gure 1 shows a comparison between these two functions for several values of the parameter q. 
Clearly, the correlation function decreases more rapidly as q increases from the standard value q = 1- At the critical 
temperature, ^ — > cx) and the correlations become 

Gc(R) « (24) 

which, due to the inequality Ac > Af ^, is smaller than G^^ . Figure 2 shows a for T > Tcriucai and some values of q 
in comparison to the q = 1 case. Since, at T > Tc, a is smaller for g = 1 we find that the correlation function for all 
temperatures is larger for the q = 1 case. 

B. Critical behavior 

In order to study the behavior at the critical temperature we need to expand the function 

2 yV^ 

[l + y-{q-l)\nzY"-'^-'' -I 

in powers of a. It is clear from the previous discussion that 63/2(2:, q) becomes the well known Bose-Einstein function 
93/2(2) in the limit q ^ 1. For these purposes we apply the same method used in Ref. p^ , which gives a power 
series of a by performing first a Mellin transformation of the function and then applying the corresponding inverse 
transform. The Mellin transform of 63/2(2, q) is given by the equation 



/ y 

63/2(^,9) = ^7 / : ,ii(„-i) dy, (25) 



^3/2(5) = / 63/2(2,9)0" Ma, 



r(5) ^ r(m/(g-l)-(3/2)-g) 

2^ v^m/(n-^\^ ' ^^^> 



_ 1)3/2+. r(m/(Q-l)) 



and the inverse transformation 



r ( \ ^ 1 ^ \- r(m/(g-l)-(3/2)-.) _ 

where the contour of integration closes on the left half plane. With use of the approximation ll^ 
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r(m/(g-l)-(3/2)-.) 



T(m/{q-l)) 



(m/(g- l))3/2+^- 



m/{q-l) {m/{q-l)Y 



(28) 



with the constants 



5=i(3/2 + s)(3/2 + ,s + l), 



(29) 



the function £'3/2(27'?) is written 



r(s)a- 



00 (g- 1)^/2+ 

+ c(g-l)^/2+^C(7/2 + s) + ...l 



{q - 1)3/2+^(3/2 + s) + 6(g - l)5/2+^C(5/2 + s) 



(30) 



The function T{s) has simple poles at s = —n with residues (— l)"/n! and the zeta function C(''^) has a simple pole 
at w = 1 with residue equal to +1. With use of the residue theorem the function G^^/^iZjq) is expressed as a power 
series of a as follows 

CM.. .) - — ^ E E (-°(^^7'))"n-../(,- l)-3/2 + ..) , r,-l/2)oV» - i(, - l)^r(-5/2,o»/. , 0(o"»,. 



m— 1 n— 



(31) 



A simple check shows that the series with integer powers of a reduces, as 5 — > 1, to the result in Ref. |1J] 
^^^q(— a)"C(3/2 — n)/n\. At lowest order in a we obtain 



{q-lf 



1 ^ r(m/(g-l)-3/2) , 



(32) 



The series in equation ( |3^ ) approximates at lowest order to the function (^(3/2) and the last term is identical to the 
standard, q — \, result. Therefore, since for a « we have £3/2(2, q) « a^^^, defining t — {T ~ Tc)/Tc wc find that 
a ^ . In addition, the correlation length ^ ^ \t^^ , as the standard case. Since, near the critical temperature, 
the t dependence of the functions G^{z^q) is the same as the case of the Bose-Einstein functions gu{z) and the 
thermodynamic functions of the two systems have the same functional form in terms of Gu and g^, the remaining 
critical exponents are also independent of q. 



III. CONCLUSIONS 



In this paper our main concern has been to study the long range behavior predicted by thermodynamical systems 
described by a factorized partition function . As pointed out in the Introduction, this factorized partition function 
has been shown to approximate well the partition function of nonextensive statistics. A calculation of the correlation 
fimctions indicate that a Bose gas for the q ^ 1 case is less correlated than for q = 1. In particular, we showed 
that the correlation length for g = 1 is larger than for q ^ 1 dX all temperatures , and the critical exponents are 

6 



independent of q. Certainly, our calculations can draw conclusions on the long range behavior of a thermodynamical 
system in the factorization approach only. On the other hand, correlations in nonextensive thermodynamics should 
be studied with use of the Tsallis' partition function. The fact that a Bose-Einstein gas obeying the studied factorized 
partition function is less correlated than for the Boltzmann-Gibbs case implies that the long range behavior expected 
in nonextensive thermodynamics is lost when the Tsallis' partition function is replaced by a factorized one. Thus, our 
results point out that the errors introduced by forcing factorization, which could be not significant for the evaluation 
of thermodynamic functions, are large enough that the long range behavior of Tsallis thermostatistics cannot be 
studied within this factorization approach. 
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FIG. 1. The function / = 277^2 (n)V3^^ ^here (n) is the average number of particles per volume and ^ is the 
correlation length at high temperatures for the cases of g' = 1,6/5,3/2 as a function of T/TcriUcai- TcriUcai refers to 
the critical temperature for each value of q. This graph shows that the correlation length decreases as the value of q 
increases. 

FIG. 2. The behavior of a = — as a function of the temperature for the cases q = 1,6/5,3/2, indicating that for 
low temperatures a Bose- Einstein gas is more strongly correlated for q = 1. 
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